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ABSTRACT 

Normandeau et. al. (1996) have proposed that W4 is a galactic chimney, the 
only chimney to-date identified in our Galaxy. Using the recent ~ 1' resolution 
IGA and DRAO CGPS galactic plane surveys we analyze the far-infrared and 
radio structure of the W4 chimney/supershell. We show W4 has a swept-up par- 
tially ionized shell of gas and dust which is powered by the OC1 352 star cluster. 
Analysis of the dust column density establishes there is dense interstellar material 
below the shell, directly showing the dense material which caused the lower shell 
expansion to stall. Due to much lower densities above the Galactic plane, the 
upper W4 shell achieved ’’breakout to form a Galactic chimney. Although the 
shell appears ionization bounded, it is very inhomogenous and an ionized halo 
provides evidence of significant Lyman continuum leakage. A large fi action of 
the OC1 352 cluster photons escape to large distances and are available to ionize 
the WIM (warm ionized medium) component ol the interstellar medium. 

Subject headings: Galaxy:strueture — HII regions — infrared: ISM: continuum 
— ISM:bubbles — ISM:individual(W4 Supershell) — ISM: structure — radio 
continuum: ISM 


1. Introduction 

The IRAS Galaxy Atlas (IGA) provides high spatial resolution processed IRAS 
images of the Galactic plane at far-infrared wavelengths (Cao et al. 1997). The IGA, in 
conjunction with the DRAO HI line/21 cm continuum (Normandeau et al. 1997) and 
FCRAO CO(l-O) line (Heyer et al. 1998) Galactic plane surveys, both with similar (~ 1') 
resolution, provide a powerful venue for studying the interstellar medium and large scale 


structure in our Galaxy. 


Figure 1 shows a panoramic view encompassing the W3, W4, and W5 HII regions 
(Westerhout 1958) in the Outer Galaxy. In the infrared W4 and W5 display pronounced 
loop structure, suggestive of wind-swept shells powered by OB clusters, while \V3 s compact 
size indicates a young HII region which is related to star formation activity in the VV3 
molecular cloud. Comparison of the 60 pmemission (left panel) with the 21 cm continuum 
emission (right panel) illustrates the close correspondence of the far-infrared to radio 
continuum emission. Previous studies show the ratio ol far-infrared to 2.7 GHz radio 
continuum flux density R is high for IIII regions (~ 1000) as compared with supernova 
remnants (< 10) (Furst et al. 1987; Haslarn k Osborne 1987). Based on this discriminator 
between thermal and nonthermal regions, most of the structure in Figure 1 is likely due to 
thermal emission. 

The W4 HII region provides a striking example of dynamically formed features in the 
ISM. Normandeau, Taylor, and Dewdney (1996) propose W4 is a superbubble powered by 
the OC1 352 (IC1805) OB cluster (/ = 135 deg, b = 1 deg) which has achieved breakout to 
form a Galactic chimney. The evidence for a chimney includes a large void in the HI gas 
above the IC1805 cluster (Normandeau et al. 1996) (indicated schematically in Fig. 1), as 
well as an accompanying Ha shell (Dennison et al. 1997). Directly below OC1 352 and 
prominent in Figure 1 is the smaller W4 loop (~ 40'), visible in infrared, radio continuum 
and Ha emission. The YV4 loop is interpreted as the lower half of the W4 supeibubble, 
which is smaller in size below OC1 352 because it has expanded into a higher ambient 
density medium. Basu, Johnstone, k Martin (1999) model the W4 region as wind-blown 
bubble in a stratified interstellar medium. 

The correspondence in Figure 1 between the 60 /miemission (left panel) and 21 cm 
continuum emission (right panel) is very strong. Notable is the large 100 pc loop of ionized 
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material which is mirrored in the 60 /imdata in the l.b x 2.o field. Marked neai the 
top is the OC1 352 OB cluster, the ionization source for the HII region. The ratio of the 
60yum/100/mi data (not shown) displays a gradient which locates the hottest dust near OC1 
352, decreasing to cooler dust temperatures in the loop. 

The stellar properties and distance to the OC1 352 cluster have been measured by- 
Massey, Johnson, k DeGioia-Eastwood (1995). A trio of massive stars (04-05) provide 
most of the ionizing radiation. The FIR and radio emission are not symmetric around the 
cluster: to the north (outside the boundaries of Fig. 2; but see Fig. 1) lies little emission 
except for a dense cometary-shaped cloud which is being photoionized and eroded by the 
OB cluster (Heyer et al. 1996b). This suggests much lower densities on average, north of 
the cluster, compared with the southern region shown in Figure 2. 

There is a halo of faint ionized gas around W4. The correspondence of 60 f.im and 21 
cm continuum continues to much lower intensities than can be seen in the stretch piovided 
in Figure 2. The data indicate there is extensive faint ionized gas centered on the HII 
region, but lying well outside the W4 loop. 

The picture supports the view that the OB cluster in a central role. Strong winds 
from the cluster are probably responsible for sweeping up the 100 pc sized loop (or shell) of 


material. 


2. Data and Analysis 


3. The Structure of the W4 Loop 
3.1. Radio versus Infrared Comparison 

Figure 2 demonstrates there is a strong correspondance between the 60 /rmand 21 cm 
continuum images. The correspondance is not automatic, as the fai-infraied emission is 
sensitive to warm (20 - 60 K) dust, whereas the 21 cm continuum traces ionized (> 3000 
K) gas. There is an overall correlation with slope R = hicml hotan ~ 300 with significant 
scatter. To understand the source of the “scatter”, we map the outlying points points 
back into image space, to look for geometrically distinct regions. We adopt the excess map 
technique of Gaustad k Van Buren (1993) and form a difference image, I 2 1 cm ^ f 60 inn i 

which serves to eliminate structure common to both images. 

The remarkable difference image (Fig. 2, middle panel) reveals striking spatial 
differences between ionized gas and warm dust. The loop structuie displays enhanced 21 
cm emission (white) lying interior to the enhanced 60/imemission (black). The IC1805 star 
cluster is the likely source of ionization^ crosses mark the positions of the thiee earliest- type 
(04-05) massive stars. The enhanced continuum emission on the inner edge suggests the 
loop is ionization bounded. 

The far-infrared emission is generally accepted to be thermal emission from dust heated 
by the UV component of the local interstellar radiation field (Terebey k Fich 1986). On 
this basis, the enhanced FIR at the outer loop edge appears due to non-ionizing UV photons 
from the OB cluster which travel into the compressed shell where the dust optical depth is 
significant. The infrared data (see section 3.2 show a corroborating gradient in the dust 
temperature across the shell. 
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The faint halo of ionized gas outside the W4 loop, seen primarily to the left and top, is 
not consistent with an ionization bounded II 1 1 region. This further suggests the wind-blown 
shell is patchy rather than homogeneous. If the shell is patchy, i.e. the ionization is density 
bounded at some but not all points, then ionizing photons can leak through and potentially 
ionize the low-density ICM out to quite large distances. 


3.2. Dust temperature and optical depth from infared data 

We express the dust opacity as k = (O.lcm^g 1 (250/rm/A #im (Hildebrand 19b3) 
where is the wavelength in microns and [3 = 0 - 2 is a standard range; we adopt 

0 = 1 for far-infrared wavelengths. The dust optical depth r = up where p is the gas 
density. The intensity due to thermal dust emission is then, for optically thin emission, 

1 = KpB u {T D ) where B„ is the Planck function and T D is the dust temperature. The dust 
temperature is derived from the 60 and 100 //mintensity ratio, where the 60 g ni data is 
beam-matched to the 100 pmresolution. The background emission from cold (~ 20 K) 
cirrus emission associated with HI clouds is subtracted from both the 60 and 100 //mi n i ages . 
The background model multiples the integrated HI brightness temperature image by an 
HI-to-FIR conversion constant. The adopted conversion constant was derived by Terebey 
& Fich (1986) for the infrared cirrus in a nearby (/ = 125 °) field. The background model 
has little spatial structure and is constant to within 10% over the field of interest, having 
a value 10.6 ± 1.0 MJy ster -1 at 60 /rm). The IRAS intensity ratio is then color-corrected 
(IRAS Catalogs and Atlases: Explanatory Supplement 1988); for example changing the 
median 60 to 100 intensity ratio in the loop from 0.52 to XX for (3 = 1 dust opacity. 
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3.3. Shell Model versus Ring Model 

Two physical models are suggested by the loop morphology of W4: a thin limb- 
brightened shell, or a physical ring, i.e. torus of material. The ring model predicts no 
material projected toward the geometric center. The thin shell model with thickness T 
has a path length 2T toward the geometric center and a contrast between center and 
limb-brightened edge depends on shell thickness. 

The spatial structure in the optical depth map implies there are large density variations 
across the W4 loop. Nevertheless the loop structure is sufficiently regular to attempt fitting 
a shell model. The spherical assumption is admittedly simplistic, as the IC1805 powering 
star cluster lies at the top edge, rather than geometric centei of the model. F igure 3 and 
4 display two models which give reasonable visual fits. Figure 3 displays the optical depth 
data in green whereas Model 1 is overlaid in red. Three slices through the data, are displayed 
in Figure 4. Two models are overlaid; the shell thickness in model 2 is one quarter that 
of model 1, and exhibits correspondingly stronger limb-brightening. Model parameters are 
given in Table 1. 


4. Discussion 

We have shown there is a strong correspondence between the FIR and radio continuum 
in the IRAS IGA and DRAO 21 cm Galactic plane surveys (Fig. 1). The emission is 
due to thermal processes (dust and ionized gas). Analysis of a subregion centered on the 
lower W4 loop supports the Normandeau et, al. (1996) W4 chimney hypothesis where a 
swept-up shell is powered by the OC1 352 OB cluster. Fits to the far-infrared data reveal 
the shell. If most of the shell is neutral, then the derived shell mass is about 10,000 M0 
with density 4.8 cm -3 , implying about 2.8 cm -3 for the ambient density. The ionized gas 
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appears interior to the dust, suggesting the shell is ionization bounded. However the shell 
is very inhoniogenous, varying by at least factors of three in density. A halo of ionized gas 
and heated dust is seen outside the shell on left side. This is direct evidence of Lyman 
continuum photon leakage. Analysis of the emission measure shows that even though the 
shell is ionization bounded in some areas, that most of the ionizing photons escape the 
shell. This supports recent proposals that 0 stars provide the ionizing photons which are 
needed to produce the extended WIM which is observed in external galaxies (Ferguson et 
al. 1996; Oey & Kennicutt 1997). 
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Table 1. Parameters of W4 Dust Shell Model 0 . 


Model 

/cen(°) 

i>cen{°) 

MtO 

Rtn(') 

R«» ( P ^ ) 

Teen 

nr/ (cm 3 ) 

M, fc (MO) 

i 

135.02 

0.42 

10 

30 

19 

0.00011 

4.8 

6700 

ii 

135.02 

0.42 

2.5 

30 

19 

0.00011 

19 

5300 


a Assumes D = 2.2 kpc, k =0.25 cm 2 g~‘ at 100 /on, mean mass per H = 1.4 


Table 2. Parameters of W4 Ionized Shell Model 0 . 


Model 

/cen(°) 

t>cen{°) 

SRC) 

Rm(') 

Pin (pc) 

n 2 e L( pc cm 6 ) 

n e (cm 3 ) 

M, fc (M0) 

i 

135.02 

0.42 

10 

30 

19 

270 

3.3 

4800 

ii 

135.02 

0.42 

2.5 

30 

19 

270 

13 

2600 


a Size and thickness of model shell is the same as Table 1. 
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Fig. 1. — Panoramic far-infrared and radio images of the W3-4-5 Outer Galaxy HII regions. 
The strong correspondence between the infrared and radio continuum emission toward the 
HII regions demonstrates a tight, coupling between warm dust (60 and 100 p m) and ionized 
gas (21 cm). Inverse color table used, so bright = black. Left panel shows IRAS 60 pm data; 
100 pm looks nearly identical. Notice both W4 and W5 exhibit loop structure, suggestive of 
limb-brightened shells. Display uses 1.0 < log I 60 (MJy ster' 1 ) < 2.5. Right panel shows 21 
cm continuum emission; l2 cm continuum looks nearly identical. The W3, W4, and W5 HII 
regions and HB3 supernova remnant are identified. VY3 is a region of active star formation. 
Solid line shows the approximate outline of the proposed W4 Galactic chimney; this model 
invokes a density gradient to explain breakout to the top (Normandeau et al. 1996). Display 
uses 0 < log I 2 j(K) < 1.2, in brightness temperature units. 

Fig. 2. — A closeup of the W4 loop region shows striking spatial differences between the 
infrared and radio continuum emission. Difference images shows ionized gas (white) lies 
interior to warm dust (black) across the shell. Positions of the three earliest O stars in the 
IC1805 OB cluster are marked. Default color table used, so bright = white. Left panel 
shows 60 pm data. Display uses 0.5 < log I 6 o(Jv arc min -2 ) < 2.5. Right panel shows 21 cm 
continuum data, -2.1 < log I 2 i(Jy arcmin" 2 ) < -1.2. Middle panel shows difference image, 
as described in text. Display uses -0.02 > I*// < 0.015. 

Fig. 3.— The optical depth and temperature across the W4 loop. Left panel shows back- 
ground subtracted 60 pm image. As the other panels demonstrate, most of the stiuctuie 
represents column density rather than temperature variation. The IG1805 OB clustei is the 
likely source of ionization and heating. Display uses 0.5 < log fi;u( M Jv ster ) < 2.5. 1 he 
r 100 optical depth map (middle panel) reveals dense relatively cool material below (6 < 0°) 
the limb-brightened shell. This dense material is predicted by the W4 chimney model, to 
explain why breakout occurs above rather than below the powering IC1805 cluster. The 



shell is inhomogenous; low column density extends from the lower-middle (135.0, 0.2) to the 
left. Display uses -5.0 < log t 100 < -3.0. Low intensity regions (I 60 < 0.5 MJy ster -1 ) are 
masked out to minimize unreliable structure due to noise (middle and right panel)(???check 
value). Right panel shows dust color temperature derived from the Ieo/Roo intensity ratio. 
Overplotted are Ti 00 = 0.00 02 (thin), 0.0005 (thick) optical depth contours. The dust tem- 
perature is nearly constant at 36 K, roughly centered on IC1805; to the left the wann dust 
extends far past the barely-visible shell edge. 1 he shell is much thickei to the lower right, 
there is a corresponding abrupt temperature decrease at the shell s inner edge. Display uses 

25 < T C (K) < 45. 

Pig 4 _ — a ) The W4 thick-shell model (red) is superimposed on the optical depth map 
(green) of the dust at 100 microns. The numbered lines indicate the positions of plotted 
slices, which show optical depth versus distance from the nominal shell center, b) The data 
(solid-line) and two models from Table 1 are shown. Models have been convolved with the 
100 /rm beam. Deviations from sphericity favor the thick-shell model, while the sharp horns 
of the line profiles favor the thin-shell model. The shell is patchy, but appears consistent 
with an inhomogenous shell rather than a thin ring, i.e. torus. 

Fig. 5.— a) The W4 thin-shell model (red) is superimposed on the 21 cm radio image (green) 
of the ionized gas. The numbered lines indicate the positions of plotted slices, b) The profiles 
of observed and predicted radio emission versus distance from the nominal shell center. The 
data (solid-line) and two models from Table 2 are shown. Models have been convolved with 
the 21 cm beam. The sharp horns of the line profiles favor the thin-shell model. The ionized 
shell is patchy, but appears consistent with a shell rather than a thin ring, i.e. torus. Figure 2 
shows that the ionized gas lies interior to the dust emission, suggesting the shell is ionization 


bounded. 
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